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' Abstract. The nonperturbative content of the QCD vacuum permits the occurrence 

of QCD-vaccum condensate contributions to penguin amplitudes. We calculate the 
dimension-4 < niqqq > and < agG^ > contributions to the effective Wilson coefficients 
^ , for penguin four-Fermi operators, and discuss how such contributions may contribute 

' to nonleptonic B decays. 
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O ■ Nonleptonic B decays are very important to our understanding of standard model 

physics. In particular, decays of the form B — »• hih2, where hi and /i2 are light 
hadrons vr, K, 7], rj', necessarily entail an operator product expansion framework 



p • with field theoretical coefficients Ci{fi): 

< hh2\H,fi\B >= ^J2V6KMC^{^^) < hih2mfi)\B > . (1) 

The standard model's heavy degrees of freedom associated with the top quark and 
^ I the have been integrated out of the effective Hamiltonian appearing in (1). 

d • The validity and utility of this description, as well as the effective truncation of 

higher dimensional operators from the expansion (1), must ultimately be tested 
empirically. For the case B — > hih2, it is standard to focus on the dimension-6 
four-Fermi operators Oi-e, where operators Oiu,ic,2m,2c are current-current operators 
and Oa-e are penguin operators as given below, with q being {s,d}, respectively 
for b {s, d} transitions[L = 1 — 75, i? = 1 + 75]: 

Olu = qal^LUaUfS^f^Lbfj , O2U = qaYLUi3U^^f,Lba , 

Oic = qal'^LcaCp^f^Lbp , O2C = qal^LcpCp-f^Lba , 
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(2) 



Methodologically, there exists a genuine issue as to how to partition radiative 
corrections between coefficients and the matrix elements in (1). We take the usual 
approach of allowing the matrix elements < hih2\0i{fi)\B > in (1) to be evaluated 
at tree order, with all higher order diagramatic perturbative corrections incorpo- 
rated in Wilson coefficients Ci{ii). 

We are particularly concerned with condensate contributions to penguin dia- 
grams, in which an intcrmcdiatc-statc quark- ant iquark pair emits an off-shell gluon 
that creates a final-state quark-antiquark pair. Such contributions occur by virtue 
of the nonperturbative content of the QCD vacuum, which necessarily entails non- 
vanishing condensate-dependent expressions for QCD-vacuum expectation values 
of normal-ordered products of fields, quantities which would ordinarily vanish if 
the vacuum were purely perturbative [1]. Thus the purely perturbative diagram 
(Fig. 1) must be augmented by diagrams sensitive to the nonperturbative vacuum 
expectation values of quark-antiquark pairs (Fig. 2) and multiple gluon fields (Fig. 
3). 

The purely perturbative Fig. 1 contributions to penguin operators are already 
known to next-to-leading (NLL) precision via naive dimensional regularization 
within an MS context [2] . We choose to incorporate such corrections entirely within 
the Wilson coefficients Cs-e, which are now denoted to be "effective" Wilson coef- 
ficients as follows: 
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AFi(^) = -4 / dxx(l - x)ln[l - zx(l - x)] 
Jo 
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An important feature of such purely perturbative corrections to penguin operators 
is their dependence not only on the internal quark mass rriq/, but also on the 
value of k"^ characterizing the final-state quark-antiquark pair. This behaviour is 
particularly important for CP violation, since AF{k'^/m'^,) acquires an imaginary 
part when k'^ > Aw?^,. Such an imaginary part has already been suggested as the 
origin of the strong phase necessary for the generation of direct CP asymmetry [2] . 

The occurrence of A;^-dcpcndcnt effective Wilson coefficients is also a feature of 
nonperturbative corrections arising from QCD vacuum condensates. Such conden- 
sate effects, which are already understood via their QCD sum-rule duality to the 
hadron spectrum, enter radiative corrections to penguin-mediated processes via 
Figs 2 and 3. To obtain the quark-condensate contribution to such processes (Fig 
2), we proceed by replacing the usual perturbative internal quark propagator {S^) 
in Fig. 1 with the full quark propagator S{p): 

5(p)=5P(p)+5^P(p). (9) 

The nonperturbative contribution to the quark propagator [S^^) is just 

iS^^{p) = j d^xe'P-'' < n\ : *(x)*(0) : |Q > , (10) 

where 1^2 > is the nonperturbative QCD vacuum [3]. The nonlocal vacuum expec- 
tation value (vev) in eq. (10) can be expanded in terms of local condensates. We 
note that the nonlocal vev of two quark fields does not contain a gluon conden- 
sate component [4] . The quark condensate projection of the nonperturbative quark 
propagator is taken from Ref . [3] : 

i5^P(p) = (27r)^(|^ + m,)F(p) , (11) 

where the Fourier transform of F{p) is expressed in terms of Bessel function 



/•^--"'^W^-^^^. (12) 



with the following mass-shell property: 

{p'-ml)F{p)^0 . (13) 

The dimension-3 quark condensate < qq > is defined to be the vacuum expectation 
value of the normal ordered local two-quark fields, i.e. 

< qq >-< n\ : ^(0)^(0) : |0 > . (14) 

Using the Feynman rule of the eq. (11), one can obtain the nonperturbative < qq > 

contribution to the effective Wilson coefficients in a straightforward manner. The 
relevant Feynman diagram is illustrated in Fig. 2 where the nonperturbative quark 



propagator S^^ is depicted by a disconnected line with two dots. Here we con- 
centrate on the loop portion of Fig. 2 which differs from the purterbative case of 
Fig. 1. Aside from the color factor, the vector current correlation function can be 
written as 



Tr[(j) - ^ + m)7^(|^ + m)'^y\F{p) 
{p — k)'^ — m? + ie 



(15) 



where the factor 2 in front is due to two possible insertions of the nonperturbative 
quark propagator in the fermion loop. By contracting p^^p'^ into H^J^^ and using 
the identity [3] 



j d^pk.pFip) = ilim-^j d^pe-'^^-^Fip) 
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where the second line is derived from eq. (12), one can show explicitly the transver- 
sality of the correlation function in eq. (15). Consequently, one finds upon con- 
tracting g^" into n<J^> and imposing the on-shell constraint (13) that 



n<'"''>;(A;') = 24(2mJ + A;') J 



F{p) 



k^ — 2k.p + ie 

The integrals appearing in eq. (17) are evaluated as follows [3]: 
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where the final line is derived by utilizing a tabulated integral [5]. Using the 
transversality of H^^^^, the result for the RHS of (15) is easily obtained: 
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This result is, in fact the < qq >-contribution to the vector-current correlation 
function first derived in ref. [6]. The corresponding gluon condensate contribution 
to the vector-current correlator, which may be derived by several different methods 
[3,6], is given by 



where 
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and where [6,7] 
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This contribution diverges sharply at 2; = 4, a divergence which is not removable 
by the operator realignment suggested in ref. [8]. 

We see that both < qq > and < > contributions to penguin operators occur 
via their known contributions to the vacuum polarization functions [eqs. (19) and 
(20)] within Figs 2 and 3. in exactly the same way that the purely perturbative 
vacuum polarization within Fig. 1 also enters penguin operators. This latter con- 
tribution is responsible for the vacuum-polarization Feynman loop integral AFi, 
as defined by (8), which appears explicitly in the penguin-operator correction AC 
[eq. (7)]. Consequently, the aggregate effect of these lowest-dimensional conden- 
sate contributions to penguin operators is obtained by retaining equations (3-7), 
but with AFi redefined so as to include < qq > and < > contributions to the 
vacuum polarization: 



AFi{z) = -4 / dxx{l - x)Ln [1 - zx{l - x)] 




(23) 
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Eq. (23) is expressed in terms of the renormalization-group invariant condensates 

(to one-loop order in whose magnitude are known from QCD sum-rule appli- 
cations. For example, if z is much larger than one, as would be the case for very 



light quarks [eq. (21)], the nonperturbative contributions to (23) arc approximately 
given by 

A z^nonperturbative _ IBtT^ < rUggq > _ 277 < a^G^ > 

Using phenomenological values < niqqq >— —flm^/A [9] and < agC^ >— 0.045 
GeV^ [10], these nonperturbative contributions arc seen to be three orders of mag- 
nitude smaller than corresponding purely perturbative contributions to AFi. 

Indeed, for B — > hih2 decays, we may anticipate that k"^ is typically between 
and ml/ 2 [11]. This range, however, suggests sensitivity to the z = 4 singularity 
mentioned above when the intermediate-state quark q' in (7) is identified with 
the charmed quark [recall that z in (23) can be identified with k'^/m'^, in (7)]. One 
cannot, therefore, discount the possibility that the gluon-condensate contribution to 
the penguin amplitude may be comparable to the purely perturbative contribution. 

This result has possible ramifications for investigations of CP asymmetry in 
hadronic B decays. The gluon condensate contribution to (7) [via (23)] develops 
an imaginary part as A;^ 4m^ from above that would dominate the strong phase 
of the effective Wilson coefficients (3-6). Interestingly, the < mqqq > contribution 
to (23), though small, has an imaginary part over the range < z < 4, suggestive 
of an additional source of CP violation stemming from light quark-antiquark-pair 
intermediate states. This latter absorptive contribution has been discussed in other 
contexts as a manifestation of the Goldstone theorem [12]. The singularity of the 
gluon-condensate contribution to penguin operator coefficients at z — 4, however, 
suggests a resonance interpretation. 

Consider, for example, the dileptonic B decays B — >■ Xsi~^i~. Such decays ac- 
quire a continuum charmed quark loop contribution, analogous to the quark-loop 
in Fig. 1, in addition to a resonance contribution associated with the cc resonances 
J/t/j, t/j', ... [13]. This resonance contribution can be modeled by the presence of a 
resonance propagator 1 / (M^ — k^ — iMF) in the long distance amplitude, M and 
k being respectively the resonance-mass and momentum-transfer to the the lepton 
pair. The presence of the resonance propagator leads to a peak at k'^ = domi- 
nating the branching ratio for this decay mode. For the S — > /11/12 case, the z — 4 
gluon-condensate singularity of AFi within AC [eq. (7)] suggests the occurrence 
of a "charmonium" resonance at k"^ = 4m^. However, the fact that this interme- 
diate state couples to a gluon, as opposed to the electromagnetic-penguin photon 
appropriate for B — > Xgi'^i^ , necessarily implies that this cc resonance be a weakly 
bound colour octet state. We note that the idea of colour-nonsinglet intermediate 
states within QCD is quite standard; such states are, of course, quarks and glu- 
ons, which are prohibited from manifesting themselves as physical final states by 
additional dynamics (confinement). What is new here is the idea of a (similarly 
confined) colour-nonsinglet bound state also contributing to a nonleptonic process 
via QCD's nonperturbative content, the QCD-vacuum condensate. 

We reiterate, however, that the qq momentum-transfer k'^ is not an observable 
in the B —>■ hih2 decay process, in contrast to the observable dilepton momentum- 



transfer in 5 — > Xgl'^l'. It should be noted, though, that Ami ^"^^ unrea- 
sonable choice for the value of A;^ characterizing B — > /ii/i2, and that the gluon- 

condensate component of the decay amplitude may well play an important role 
both in enhancing such decays and in providing a mechanism for CP violation. 
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